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Abstract. We investigate the regularity of bounded weak solutions of scalar 

conservation laws with uniformly convex flux in space dimension one, satisfying 

an entropy condition with entropy production term that is a signed Radon 

1/3 3 

measure. We prove that all such solutions belong to the Besov space BJ^ j'^^. 
Since C. DeLellis and M. Westdickenberg [Ann. Inst. H. Poincare Anal. Non 
Lineaire 20 (2003), 1075-1085] have proved the existence of such solutions that 
do not belong to B^'f^^ if either s > 1/ max(p, 3) or s = 1/3 and 1 < g < p < 3 
or s = l/p with p > 3 and q < oo, this regularizing effect is optimal. The 
proof is based on the kinetic formulation of scalar conservation laws and on 
an interaction estimate in physical space. 



1. Introduction 

Consider the Cauchy problem for the free transport equation with unknown 
u = u{t, a;) G R 



(1) 



dtU + cdxU — 0, a;£R,<>0, 



where c G R is a constant. It is well known that 

u{t,x) = M™(a; - ct) , 

so that u{t, •) has exactly the same level of regularity as w™. 

If the speed of propagation c depends on the unknown u, the situation is com- 
pletely different. Consider the scalar conservation law with unknown u = u{t, x) £ 
R and flux a : R — > R of class 



(2) 



dtu + dxa{u) = 0, a;eR,t>0, 



If u is of class C^, the scalar conservation law (2) is equivalent to the free transport 
equation (1) with c = a'{u). But even if u"' G C^(R), it may happen that the 
solution u of (2) loses the regularity in finite time. (This was already known 
to Riemann: see [34].) More precisely, if a is convex, and if m*" is decreasing on 
an open interval, there exists T G R!^ such that (2) cannot have a solution 
defined on (0, T') x R for all T' > T. However, (2) has global weak solutions whose 
restrictions to (T, oo) x R contains jump discontinuities. 
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If a"{v) > a > for each u £ R, for each m'" e L^(R), there exists a unique weak 
solution u G L°°(R+; L^(R)) of (2) satisfying in addition the differential inequality, 
referred to as the entropy condition 

(3) dt-niu) + < 

for each convex function rj, referred to as the entropy, where 

(4) q{v) = j T]' {'w)a' {'w)d'w , 

referred to as entropy flux associated to t}. 

This solution satisfies u{t^ •) £ BVioc H _L°°(R) for each t > — see formula 
(4.9) in chapter 4 of [26] for the L°° bound, and chapter 16 §A in [37] for the 
BVioc bound. Thus, if u*" S (7^(R), since u{t,-) € BVioc(R) niay contain jump 
discontinuities, the effect of the nonlinearity a is a loss of regularity in the solution. 
Yet, if w*" € i^(R-), the fact that u{t, •) e BVioc(R) for each t > can be viewed 
as a (limited) regularizing effect. 

The purpose of the present paper is to study the optimal regularizing effect for 
weak solutions of (2) satisfying the weaker entropy condition 

(5) dtr]{u) + dxq{u) = -/i 

for each entropy-entropy flux pair {rj,q) as above, i.e. satisfying (4), where is a 
signed Radon measure on R^ x R (instead of a positive measure). Such solutions 
may contain jump discontinuities that would dissipate instead of create entropy and 
therefore would be considered as unphysical by analogy with gas dynamics. Yet, 
such solutions are relevant for other physical applications, such as micromagnetism: 
see for instance [35] and the references therein. 

Another motivation for considering the entropy condition (5) with entropy pro- 
duction n of indefinite sign can be found in the work of Hwang and Tzavaras [21]. 
In this work, the authors show that two different strategies for approximating so- 
lutions of scalar conservation laws, the relaxation approximation d la Jin-Xin [24] 
and the diffusion-dispersion approximation d la Schonbek [36] , may lead to kinetic 
formulations involving entropy production measures that may in general fail to be 
positive. 

The regularizing effect for this type of solutions of (2) has been studied so far by 
using a kinetic formulation of the scalar conservation law: see [27] for the original 
contribution, [22] for an improved regularity result, and [32] for a detailed presen- 
tation of kinetic formulations and their properties. The tool for establishing the 
regularizing effect for kinetic formulations is a class of results known as velocity av- 
eraging, introduced independently in [1] and [19], with subsequent generalizations 
and improvements described for instance in chapter 1 of [7] — see also the list of 
references given in section 3. 

The best result obtained by this method is that u G W^^^ x R) for all s < 1 /3 
and 1 < r < 3/2 (sec [22]) in the earlier result [27], the integrability exponent 
was restricted to r < 5/3. On the other hand, it is possible to construct such 
solutions that do not belong to any Besov space "better than" B^^i^^^RZ^ x R) — 
see section 4 for a more precise statement of this optimality result and [11] for the 
proof. 

In the present paper, we prove that solutions of (2) with entropy production 
that is a signed Radon measure for each convex entropy does indeed satisfy a 
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(local) Boo ' estimate provided that the flux function is uniformly convex. With 
the result in [11], this shows that the optimal regularity space for such solutions is 
indeed B^J^^iiRX x R). 

Theorem 1.1. Assume that a £ C^(R) satisfies a" > ao > on R, and let 
It*" G L°°(R). Any bounded weak solution of (2) satisfying (5) for each convex 
entropy rj, with an entropy production fi that is a signed Radon measure on R+ x R, 
belongs to Bl^^^^iiRl x R). 

Our proof of this result relies on a method completely different from velocity 
averaging. At variance with velocity averaging, this method does not use any 
argument from harmonic analysis (Fourier transform, Littlewood-Paley decompo- 
sitions...), but is based on an "interaction identity" presented in section 2. After 
recalling an earlier, weaker regularizing effect obtained in [15, 16] for solutions of (2) 
subject to a weaker entropy condition — Theorem 5.1 in section 4 — the optimal 
regularizing effect is stated as Theorem 4.1 in that same section. As a warm-up, 
we use our method based on the interaction identity to establish a new velocity 
averaging result — Theorem 3.1 in section 3. Unlike in all velocity averaging the- 
orems known to this date, the regularity of velocity averages of the solution of the 
transport equation obtained in this result is independent of the regularity of the 
source term in the velocity variable. The tradeoff is that the kinetic solution must 
satisfy a seemingly unnatural monotonicity condition — see condition (11) below 
— which however turns out to be satisfied precisely in the context of the kinetic 
formulation of scalar conservation laws. 

This new method is reminiscent of some tools from compensated compactness 
[28, 39] — see section 2 below and the comments in section 2 of [15] . A shortcoming 
of this method is that, unlike velocity averaging, it seems confined to the case of one 
space dimension, at least at the time of this writing. In the case of space dimension 
higher than one, some regularizing effect in the time variable of special classes of 
nonlinear fiuxes (that are in particular power-like at infinity) has been established 
in [30] — an earlier result of same type in the case of homogeneous nonlinearities 
can be found in [3]. 

2. The Interaction Identity 

Consider the system of partial differential relations 

f dtA + d^B = C, 
{dtD + d^E = F, 

where A, B, C, D, E, F are real-valued functions of t > and a; G R. We henceforth 
assume that the functions A, B, C, D, E, F have compact support in R^ x R, and 
are extended by to R x R. 
The quantity 




A{t, x)D{t, y)dxdv 



has been introduced by S.R.S. Varadhan — see Lemma 22.1 in [41] — in the context 
of discrete velocity models in kinetic theory, and used by several authors since then 
(J.-M. Bony [5], C. Cercignani [8], S.-Y. Ha [20]). It is also reminiscent of Glimm's 
interaction functional used in the theory of hyperbolic systems of conservation laws 
in space dimension 1: see [14]. 
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Differentiating under the integral sign and using both equations in (6), one ob- 
tains the following interaction identity 



/ / iC{t, x)D{t, y) + A{t, x)F{t, y))dxdy 
- 11 {d^B{t,x)D{t,y)+A{t,x)dyE{t,y))dxdy 

J Jx<y 

= jj iC{t,x)D{t,y) + A{t,x)F{t,y))dxdy 

J J x<y 

+ [ {AE-DB){t,z)dz. 

Integrating further in the time variable and observing that 7i is compactly sup- 
ported in t e R^, we arrive at the identity 



(7) 



/ / (AE - DB){t, z)dzdt = - [[ C{t, x) ( ( D{t, y)dy) dxdt 

- [[ F{t,y)( r A{t,x)dx]dydt. 
J Jnxn \J-oo / 



Exchanging the roles of the time and space variables in the computation above, 
one can consider the quantity 



l2{t) := B{s,z)E{t,z)dsdt 

J J 8<t 



' S<t 

instead of 7i . Proceeding as above, one sees that 



^{z)= [[ {C{s,z)E{t,z) + B{s,z)F{t,z))dsdt 
dz JJs<t 



- {dsA{s, z)E{t, z) + B{s, z)dtD{t, z))dsdt 

J Js<t 

= 11 (C{s, z)E{t, z) + B{s, z)F{t, z))dsdt 

J J s<t 

- I {AE - DB){t,z)dt. 

Integrating further in z and observing that I2 is compactly supported in 2; G R, we 
obtain 



/ / {AE - DB){t, z)dt =11 C{s, z)( [ E{t, z)dt) dzds 
JJuxii JJrxii \Js J 

[ [ F(t, z)( [ B(s, z)ds] dzdt . 



(8) i-i. _ ^ 

+ 



In his proof of the interaction identity on p. 182 of his book [41], L. Tartar ob- 
serves that the structure of this identity is reminiscent of compensated compactness 
[28, 39]. Indeed, introducing the vector fields 

U{t,x,y) ■.= {A{t,x),B{t,x),Q), V{t,x,y) = {E{t,x),-B{t,x),0) 

we see that the system (6) takes the form 

div [/ = C , 
curly = (0,0, -F). 
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While the left hand side of (7) involves the inner product U ■ V, the right hand 
side involves integrands in the form of binary products where one of the terms is 
integrated — and therefore gains one order of regularity — in the space variable. 
The same is true of (8), by which one can hope to gain one order of regularity in 
the time variable. 

In view of this observation, the idea of using the interaction identities (7)-(8) for 
the purpose of establishing regularization result appears fairly natural, and will be 
used systematically in the sequel. 

3. Velocity Averaging in Physical Space 
Let / = f{t, X, v) satisfy 

(9) {dt + a'{v)d,)f = dZm 

where a is a smooth function while m is a bounded, signed Radon measure on 
R+ X R X R, and 7 e N. 

Transport equations of this type naturally appear in the kinetic formulation 
of hyperbolic systems of conservation laws: see for instance [32]. A fundamental 
question in the context of kinetic models is to investigate the local regularity of 
moments in the velocity variable v of the function /. Systematic investigations on 
this class type of questions began with our work with R. Sentis [19] (see also the 
independent study by V. Agoshkov [1]) and in a scries of subsequent contributions 
by several authors where more and more general classes of functions / and right 
hand sides are considered: see in particular [17, 13, 6, 12, 33, 18, 38, 22, 23, 4, 2]. 

All these works use at some point tools from harmonic analysis: Fourier trans- 
form, Hardy-Littlewood decomposition, Radon transform. Moreover, in all these 
results, the regularity of moments in v of the function / depends on 7. 

In this section, we give an example of velocity averaging result where the reg- 
ularity of the moments in i; of / is independent of 7, at the expense of an extra 
assumption on the v dependence in /. Also, the proof of this result is based on the 
interaction identities (7)-(8) and uses only elementary techniques in physical space. 

Theorem 3.1. Let a € C^(R); assume that there exists /3 > 1 such that, for each 

M > 0, there exists um > for which 

(10) a'{v) - a'(w) > aM{v - wY , -M <w<v<M. 

Let 7 e N and let m be a signed Radon measure on H+ x R x R. 

Assume that f G i°°(R+ x R x R) satisfies (9) and that, for each y G R and 
each s > 

{f{t + s,x + y,v)- f{t, X, v)){f{t + s,x + y,w)- f{t, x, w)) > 

for a.e. {t,x,v,w) € R+ x R x R x R. 
Then, for each ijj G C^(R), one has 

f^l^{v)dv G i?:,^,„,(R; x R) wtth s = . 

We recall that, in the context of velocity averaging, the gain of regularity ob- 
tained on averages of the form 

/ f{t,x,Od^ 

J|{|<M 
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of solutions of the transport equation 
involves 

sup |{^ e s.t. ICI < M and \w + V{^) ■ A;| < e}| , 

where, for each A c R^, the Lebesgue measure of A is designated by \A\ — see 
condition (2.1) in [17]. Equivalently, whenever k ^0, 

\{^ e R^ s.t. Id < M and |w + F(0 • k\ < e}\ 



1^ G R^ s.t. 



Id < M and F(0 • 1^ e 



OJ + € € — U 



\k\ ' |fc| 

On the other hand, (10) implies that a is strictly convex and is equivalent to the 
condition 



\{v e R s.t. \v\ < M and a'{v) G [A,B]}\ = a'-\B) - a'-'^{A) < 



B-A 

aM 



1//3 



for each A, B such that a'{-M) < A< B < a'{M). Therefore, (10) is a condition 
of the same type as the classical condition used in velocity averaging. 

A typical sufficient condition under which a satisfies (10) is as follows. Assume 
that a G C^"(R) is convex and and that, for some n G N* and z G (— M, M), one 
has 

a"{z) = ... = 0,(2"-!) (z) = , and a^^") (^) > . 
(An example of this situation is the case of 

a(.):=i^.-, 

with z = 0.) By continuity of a^^"^ there exists p > such that [z — p, z + p] C 
(-M, M) and A > such that 

a(2n) (i) > A > for each t E [z - p,z + p] . 

By Taylor's formula, whenever z — p<w<v<z + p, one has 

a'{v) - a'iw) = £ j^^^^iiv - t)^"-'lt>z + it- wr--H,^M'''\t)dt 
-((t;-^)2"-i + (^-M;)2"-i) 
;max((u-^)2"-i,(2-w;)2"-i) 

22"-i(2n-l)!^ ' 
On the other hand, whenever — Af <w<z — p<z + p<v< M, one has 
a'{v) — a'{w) = a'iv) — a'{z + p) + a! [z + p) — a'{z — p) + a'{z — p) — a'{w) 
> a'{z + p) — a'{z — p) 
A 



> 



> 



> 



A 

(2n- 1)!' 
A 

(2n-l)! 
A 



> 



> 



22"-i(2n-l)! 
A 

22"-i(2n- 1) 



(2p)^ 



\ 2n— 1 



OPTIMAL REGULARITY FOR CONSERVATION LAWS 



7 



where the first inequahty foUows from the convexity of a, while the second is a 
consequence of the previous estimate in the case v = z + p and w = z — p. Therefore 
a satisfies (10) with /3 = 2n — 1 and aM = 2^^-^(2n-i)\ (m) " • This situation is 
essentially the same as the one described in assumption (H) before Theorem 2.2 in 
[15]. 

The assumption (11) may seem somewhat unnatural; however, as wc shall see 
below, it is relevant in the context of conservation laws. A typical example of 
function / satisfying (11) is as follows. 

Let p £ L°°(R+ X R), and let T/K : R ^ R be a nondecreasing or nonincreasing 
function; then the function / defined by the formula 

f{t,x,v)=W{p{t,x)-v) 

satisfies (11). Indeed, for each t,s > and each x,y £ li such that p{t,x) and 
p{s, y) are defined 



/.p(t,x) 

f{t,x,v) - f{s,y,v) = W'{u-v)du 

J o(s.v) 



fp(t,x) 
' p{s,y) 

SO that, assuming without loss of generality that W is nondecreasing, we see that 
p{t, x) > p{s, y) =)■ fit, X, v) > f{s, y, v) for each v e R . 

Equivalently 

f{t, X, v) > f{s, y, v) for some e R p{t, x) > p{s, y) 

f{t, X, w) > f{s, y, w) for all w e R . 

Therefore 

{f{t, X, v) - f{s, y, v)){f{t, X, w) - f{s, y, w)) > for each v,w eR 

for all t,s £ R+ and all x,y G R such that p{t, x) and p{s, y) are defined, which 
means that assumption (11) is satisfied in this example. 

We do not claim that the B^^^^ regularity obtained in Theorem 3.1 is optimal. 
Since Theorem 3.1 is not the main result in the present paper, but rather an illus- 
tration of how to use the interaction identities (7)- (8) for the purpose of obtaining 
a velocity averaging theorem where the regularity index s = is independent 

of the number 7 of derivatives in v in the source term, we have left this question 
aside. 



Proof. For /i G R, define the operators 2?f and T>^ by the formulas 
(12) 



V^(j){t, x) := (l){t + h,x)- (j){t, x) , 



V^(j>{t, x) := (l){t, x + h)- (j){t, x) . 
Pick X G ^^(R^;. X R). One seeks to estimate 



jj X{t,xf V>i j f{t,x,v)i,{v)dv 



2 

dxdt . 
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This quantity is decomposed as 

JJ X{t,xf J f{t,x,v)ij{v)dv dxdt 

X(t, xfVlf{t, v)V'^f{t, X, w)tp{v)tp{w)dvdwdxdt 

= jj jj {1 - Mw - v))x{t,xfV^f{t,x,v)V^f{t,x,w)i>{v)i>{w)dvdwdxdt 

e{w - v)x{t, xfV^f{t, X, v)V^f{t, X, w)il){v)il){w)dvdwdxdt 

= Ji + J2 , 

where (t>e{z) = ^{z/e) and $ e C°°(R) satisfies 

< ^{z) < 1 , ^{z) = 1 if 1^1 > 2 and ^{z) = if 1^1 < 1 . 
By assumption (11), the integral Ji is estimated as follows: 

\Ji\< // xit,xfV^f{t,x,v)V^f{t,x,w)ipiv)ip{w)dvdwdxdt 

JJ JJ\v-w\<2e 

<4||x'lUi||/||ioo // i;{v)i;{w)dvdw 

J J\v-w\<2e 

<Co€, 

where 

Co:= 16||x'l|Li||/||ioo||V||Li||V'l|L~. 
As for the integral J2, using again assumption (11) shows that 

IJ2I < C,J3 

where, by (10), 

0e(w-v) 1 

sup , < 



_v<w<v<v {v - w){a'{v) - a'{w)) ayei+^ ' 
assuming without loss of generality that supp(V') C [—V, V], and 

J3 = jjjj {v-w){a'{v)-a'{w))x{t,xyV^f{t,x,v)V^f{t,x,w)^{v)^{'w)dvdwdxdt. 

The integral J3 is now estimated by the interaction identity (7). 
Set 

A{t, X, v) := x{t, x)V'^f{t, X, v) , D{t, x, w) := xit, x)V'^f{t, x, w) . 

Since / = f{t,x,v) satisfies (9), the functions A and D defined above satisfy (6) 
with 

B{t, X, v) := a'{v)x{x)V^J{t, x, v) , E{t, x, w) := a'{w)x{x)V^f{t, x, w) , 
and 

J C{t,x,v) ■.= x{t,x)d2V'im{t,x,v) +X{t,x,v)V^f{t,x,v), 
\ F{t, X, w) := x{t, x)dZV^m{t, x, w) + X{t, x, w)V^f{t, x, w) , 
with the notation 

X{t, X, v) = {dtx + a'{v)dxX){t, x) . 
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(We have abused the notation m(t,x,v) as if the signed Radon measure m was a 
function.) 

At this point, we apply the interaction identity (7). After multiplying both sides 
of this identity by {v — w)'ip{v)'ip{w) and integrating in v, w, one obtains 



(14) 



J3 = — J J {v — w)il}{v)'^{w) J J C{t,x,v) (^j D{t,y,w)dy^ dxdtdvdw 

~ II ~ II -^(^'y'"*^) (^J A{t,x,v)docJ dydtdvdw . 

Each integral in the right hand side of J3 involves two different kinds of terms. One 
is 

(15) 

J31 = J J {v-w)tjj{v)tjj{w) J j X{t,x,v)V^f{t,x,v) (^J D{t,y,w)dyj dxdtdvdw, 

the other being 
(16) 

J32 = JJ {v-w)4>{v)'tl>{w) J J x{t,x)d^'D^m{t,x,v) (^j D{t,y,w)dy^ dxdtdvdw. 
In both expressions, the inner integral is put in the form 

POO poo poo 

/ D{t,y,w)dy= V'^{x{t,y)f{t,y,w))dy- f{t,y + h,w)V'^x{t,y)dy 

J X J X J X 

px-\-h poo / ph \ 

= -j Xit,y)f{t,y,w)dy - J f{t,y + h,w)ij d^xit,y + z)dzj dy , 
so that 



(17) 
Thus 



D{t,y,w)dy 



<ll/llL~(llxlU~ + l|a.xlUO|/^l- 
IJ31I <Ci|/i| 



with 

Ci = 2||/||ioo(||x||Loo + \\^,x\\L^){\\vXi;U^\mL^ + ||t;V|UO • 

In J32, we first integrate by parts and bring the v derivatives to bear on the 
weight {v — w)ip{v): 

Jz2 = {-ir jj d:^{{v-w)^{vmw) jj x{t,x)V^,m{t,x,v) 

X D{t,y,w)dy^ dxdtdvdw. 

Assuming without loss of generality that supp(x) C [0, T] x [—R, R] and recalling 
that supp('^) C [— V, 1^] while \h\ < 1, one obtains 

IJ32I <C2|/l| 

with 

C2 = 2\\f\\L^{\\d^v^)\\L^ ii^iUi + iia-'v^iiL- Ikv^lUO 

pT pR+1 pV 

x(llxlk- + l|S^xlUOIIxllL~ / / \m\. 

Jo J-R-1 J-V 
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In conclusion 
Therefore 



IJsl <2(Ci+C2)|/i|. 



V'i / f{t,x,vmv)dv 



dxdt<Coe+^{C^+C2)^ 



and choosing e = we find that 



(18) X{t,xf 



dxdt < (Co + + C2)) \h\'/^^+^^ . 



As for the time regularity, it is obtained similarly, exchanging the roles of the 
variables t and x. We briefly sketch the argument below. One seeks to estimate, 
for each h > 0, the quantity 

2 



xf 



2?f / f{t,x,v)ij{v)dv 



dxdt 



that is decomposed as 



JJ X{t,xf -DlJ fit,x,v)^jj{v)di 



dxdt 



+ 



X{t, xfV^f{t, X, v)V^f{t, X, w)ip{v)ip{w)dvdwdxdt 



(1 -(peiw- v))x{t, xfV^fit, X, v)V^f{t, X, w)ip{v)ip{w)dvdwdxdt 

e{w - v)x{t, xfV^f{t, X, v)V^f{t, X, w)il){v)il){w)dvdwdxdt 

= Kr+K2. 



The term Ki is obviously estimated exactly as J\: 

(19) lifil <4||x'||Li||/||ioo / / i,{v)i,{w)dvdwdxdt < Coe. 

J J\v-w\<2e 



As in the case of J2, one has 



\K2\ < c,Ks 



where 

Ks = jjjj {v-'w){a' {v)-a' {'w))x{t,xfV^ f{t,x,v)V^ f{t,x,w)il>{v)il>{w)dvdwdxdt. 

This last integral is estimated by using the interaction identity (8), with a slightly 
different definition of A, B, C, D, E, F: 

A{t, X, v) := x{t, x)V'lf{t, X, v) , D{t, x, w) := x{t, x)V'^f{t, x, w) , 

while 

B{t, X, v) := a'{v)x{x)V^f{t, x, v) , E{t, x, w) := a' {w)x{x)V'^ f{t, x, w) 

and 

r C{t, X, v) = x{t, x)dZV^m{t, X, v) + X{t, x, v)V^f{t, x, v) , 
1 F{t, X, w) = x{t, x)dZV^m{t, X, w) + X{t, x, w)V^f{t, x, w) . 
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Apply the interaction identity (8) after multiplying both sides of this identity by 
{v — w)il){v)il){w) and integrating in v, w: 



(20) 



Ks = — JJ {v — w)tp{v)ijj{w) JJ C{s,x,v)^J E{t,x,w)dt^ dxdsdvdw 

— JJ {v — w)iIj{v)iIj{w) JJ F{t,x,w) ^J B{s,x,v)ds^ dxdtdvdw . 

As in the case of J3, the right hand side of K3 involves two different kinds of terms: 
(21) 

K31 = J J {v-iu)ip{v)ip{w) J J X{t,x,v)T>^f{s,x,v) (^J E(t,x,w)ds^ dxdsdvdw 

and 
(22) 

K32 = JJ {v^)'4>{v)4'{w) JJ x{t>x)d^V'^m{s,x,v) ^J E{t,x,w)dt^ dxdsdvdw . 
The inner integral is put in the form 

/■OO POO 

/ E{t,x,w)dt= V^{x{t,x)a'{w)f{t,x,w))dy 

J S •/ S 

a'{w)f{t + h, X, w)V^x{t, x)dt 



- / x{t,x)a'{w)f{t,x,w)dt 

J s 

POO 

- / a'{ 

J s 



(L 



{w)f{t + h,x,w) { / dtx{t + r,x)dT\ dt, 



so that 

(23) 
Thus 

with 



E{t, X, w)dt 



<W{w)\\\f\\L^{\\x\\L^ + \\d^x\\L^)h. 
li^isl <C3/l, 



C3 = 2ii/iiioo(iixikoo + \\dtx\w){\\vxnLAWnL^ + m\\LA\va'nL^). 

Next 

Ks2 = {-ir JJd^{{v-w)i,{v))i,{w) J J x{t,x)V^m{s,x,v) 

J E{t,x, w)dt\ dxdtdvdw . 



Assuming that supp(x) C [0,r] x [—R,R] while supp(V') C [— V^, V^] without loss of 
generality, and that < h <1, one obtains 

|i^32| < C^h 

with 

Ca = 2\\f\\L^{\\d^vi;)\\L^ llaVIUi + Wd-^i^h^ |ka>|UO 

x(iixik~+iiaaiuoiixiiL~ r^' r r h. 

Jo J-rJ-v 
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In conclusion 
so that 

X{t,x) 

and choosing e = we find that 



|/^3|<2(C3 + C4)|/l|, 



dxdt<Coe + ^{C3 + C4)^ 



dxdt<(Co + l{C^ + C4))h}/^^+^y 



(24) jjx{t, 

Putting together (24) and (18), we conclude that 

j mv)dv e b:;%,(r; X R) with s = ^ 



2/3- 



□ 



4. Optimal Regularizing Effect for Scalar Conservation Laws in 

Space Dimension One 

Consider the scalar conservation law 

dtu + dxa{u) = 0, X gR, t > , 



(25) 



Assume that u is a weak solution whose entropy production rate is a signed Radon 
measure. Specifically, we mean that, for each convex entropy r] € C^(R), one has 



T]" {v)dm{-, •, v) 



R 



where m is a signed Radon measure on R+ x R x R (with compact support in the 
variable v) and the entropy flux q is deflned by (4) . 

Equivalently (see for instance §6.7 in [9]), u satisfies the following kinetic formu- 
lation of (25) 



(26) 

where 
(27) 
while 
(28) 



dtf + o!{v)dxf = dytn , x,v €11, t > , 



f\ = P' 
■I lt=o •' 



f{t,x,v) :-- 



r{x,v) := 



+ l[0,«(*,a;)](l^) 
-'^[u{t,x)fi\{v) 



+ l[0,u»"(2;)](^) 
-l[ui"(2:),0](^') 



if u{t,x) > 0, 
if u{t, x) <0, 



if u'"{x) > 0, 
if u'"(a;) < 0. 



Theorem 4.1. Let a € C^(R) satisfy (10), let m he a signed Radon measure on 
R+ X R X R. Let u'" e L°°(R) and let u e L°°([0,T] x R) satisfy (26). Then 

w e ^^/^^''.(r; X R) with p = 2 + p. 
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More precisely, for all e > and all ^ G [— e, e], one has 



<2|^|(||x||L== + ||a,xlU0 2[/||X|Ui + ||x||l~ / / / d\m\\ 

\ Jo J-R-eJ-U J 



and, for allO < t < e 

^ ' ' x{t,xf\u{t + T,x)-u{t,x) f+^dxdt 



(/3 + !)(/? + 2) Jo Jr 
< 2r(||x|koo + IIS.xIIlO l\\a'\\L^-u,u)\\X\\L^ + \\x\\l^ [^11 

\ Jo J-rJ-U 



where U := while X{t,x,v) := {dt + a'{v)dx)x{tjx) cmd T,R > are such 

thatswpp{x) C [0,T] X [-R,R]. 

The following statement is an obvious consequence of this theorem. 

Corollary 4.2. Under the same assumptions as in Theorem 4-1, one has also 

u e B^^;^^(R+ X R) for each p > 2 + /3 . 

Proof Pick X e C~(R!;_ x R), and extend xu by to R^. Then, for eachp > 2 + /3 
and T, ^ e R, one has 



// 



\xu{t + T,x + ^) — xu{t, x)\^dxdt 

< {2\\xu\\L^r-'-^ II \xu{t + T,x^£,)-Xu{t,x)\'+^dxdt<C{\T\ + \£,\) 



since xu € x B^'^^'^+'^IR^). □ 

Before giving the proof of this theorem, a few remarks arc in order. 

First, the exponent P in (10) can be viewed as a measure of the nonlinearity in 
(25). Indeed, if a is linear, a' is a constant so that (10) holds (for |t« — 1;| < 1) with 
/3 = +00. and there is no regularizing effect. In other words, the regularizing effect 
predicted in Theorem 4.1 is a consequence of the nonlinearity. 

In the case where a satisfies 

(29) a e C^(R) with a"{v) > a > , veK 

P. Lax [25] and O. Oleinik [29] have proved that, for each initial data S L^(R), 
the Cauchy problem (25) has a unique entropy solution, i.e. a weak solution satis- 
fying 

dfiiiu) + d^qiu) < 

for each convex rj G C"'^(R) and q defined by (4), and that this solution u satisfies 
u{t, •) G L°°(R) for each t > {) and together with the one-sided bound 

dxuit, x) < — , t> Q , a;GR. 
at 

A consequence of this one-sided bound is that u G BVioci^+ x R). 

In the case where a{v) = ^v"^, for each p G [1, oo] and each a > 1/ max(3,j3), C. 
DeLellis and M. Westdickenberg [11] have proved the existence of a weak solutions 
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u of (25) satisfying the entropy relation (42) for each convex ij 6 C^(R) with q 
defined as in (4) and an entropy production rate jj, that is a signed Radon measure 
on R+ X R, such that 

Of course, the condition (29) implies that (10) is verified with aM ~ ct for each 
M > and (3 = 1. In that case, Theorem 4.1 predicts that weak solutions u of (25) 
satisfying the entropy relation (42) for each convex r] € C"^ (R) with q defined as in 
(4) and an entropy production rate /i that is a signed Radon measure on R+ x R 
belong to the Besov space 

Bliloc(^*+ X R-)- According to the result of C. DeLellis 
and M. Westdickenberg [11], this regularity is optimal. 

A key argument in the proof of Theorem 4.1 is the following inequality. 

Lemma 4.3. Assume that a € C^(R) satisfies assumption (10). For all u G R, 
define 

^-^^)-\ -l[„,o](«) ifu<0. 

and 

A{u,u) := jj lR^{v-w)ia'iv)-a'{w))iMu{v)-Muiv)){Muiw)-Mu{w))dvdw . 

Then, for each V > and u,u G [—V, V], one has 

A(u u) > "^^^ \u - zZp+^ 

The proof of this inequality is deferred until after the proof of Theorem 4.1. 

Proof of Theorem 4.1. Pick x € C^{R+ x R), with support in [0,T] x [-R,R], 
and let U = \\u\\l'^(^[o^t]xii)- 

We first establish the regularity in the space variable x. As in the proof of 
Theorem 3.1, set 

A{t, X, v) := xit, x)Vlf(t, X, v) , D{t, x, w) := x(i, x)Vlf{t, x, w) . 

Since / = f{t,x,v) satisfies (9), the function A and D defined above satisfy (6) 
with 

B{t, X, v) := a'{v)x{t, x)V^fit, x, v) , E{t, x, w) := a'{w)x{t, x)V^fit, x, w) 
and 

J C{t,x,v) := xit,x)dvV'^m{t,x,v) +X{t,x,v)V'^f(t,x,v), 
[ F{t, X, w) := x{t, x)d^V^m{t, x, w) + X{t, x, w)V^f{t, x, w) , 

with the notation 

X{t, X, v) := dtx{t, x) + a'{v)dxX{t-, x) ■ 
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Multiplying each side of the interaction identity (7) by 1r_^ (v—w) and integrating 
both sides of the resulting equality in the variables v and w, 
(30) 

Q:= / / lR+(t^ — / / {A{t,x,v)E{t,x,w) — B(t,x,v)D{t,x,w))dxdtdvd'w 

= — // iR+iv — w) // C{t,x,v)i / D{t,y,w)dy\ dxdtdvdw 

~ 11 1r+(^^ — w) // F{f,y,w){ / A{t,x,v)dx\ dydtdvdw . 
J Jnxn J Jnxn \J-<x J 

The right hand side of this identity involves essentially terms of two different kinds: 

(31) 

Qi = 1 1 ln+{v-w) 1 1 X{t, X, v)V^f{t, x,v)( [ D{t, y, w)dy) dxdtdvdw , 
J J J Juxr \Jx J 

the other being 

(32) 

Qi = JJ In+iv-w) JJ x{t,x)dv'D^m{t,x,v) D{t,y,w)dy^ dxdtdvdw . 

Proceeding as in the proof of Theorem 3.1, we put the inner integral in the form 
/ D{t,y,w)dy= V^{x{t,y)f{t,y,w))dy- f{t,y + h,w)V'^x{t,y)dy 

Jx J y 'fx 

f-X + h /'OO / nh \ 

= -J X{t,y)f{t,y,w)dy - J f{t,y + h,w)lj^ d^x{t,y + z)dzj dy , 



so that 
(33) 
Thus 

with 



D{t,y,w)dy 



< (llxlk- + \\dxX\\L^)\h\l[-u,u]{w) ■ 
\Qi\<K,\h\, 



(llxIk^+ll^xXlUO JJ \\X{;;v)\\Lil[_u,U]{^)M-U,U]Hdvdw 

<2C7(||x|k~ + ||9xXlUi)||^IUi =:^i. 
On the other hand, for each h € [— e, e] 

Q2 = J J 'i-R+{v — w) J J x{t-,x)dvV^m{t,x,v) (^J D{t,y,w)dyJ dxdtdvdw 

= — // x(^) 2;)2?^w(t, a;, w) I / D{t,y,w)dy \ dxdtdw , 
JnJJuxn \Jx J 

so that 

IQ2I <i^2|/l| 

with 

/^2:=2(||x||l~ + ||S.xI|lOIIxI|l- / / \m\. 

Jo J-R-e J-U 

In conclusion, one has, by the same argument as in the proof of Theorem 3.1 
(34) |g| <2(i^i + /^2)|/i|. 
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On the other hand, according to Lemma 4.3 

i-T 



Q = I I x{t,xfA{u{t,x),u{t,x + h))dxdt 
(35) ,r , 

Putting together (35) and (34) shows that 

lOC V +' QC.lOC V // 

Now for the time regularity. FoUowing the proof of Theorem 3.1, pick < h < e 
and set 

A{t, X, v) := x{t, x)V^f{t, X, v) , D{t, x, w) := x{t, x)V^f{t, x, w) . 

Since / = f{t,x,v) satisfies (9), the function A and D defined above satisfy (6) 
with 

B{t, X, v) := a'{v)x{t, x)V'lf{t, x, v) , E{t, x, w) := a'{w)x{t, x)V'^{t, x, w) , 
and 

( C{t, X, v) := x{t, x)d^V^m{t, x, v) + X{t, x, v)V^f{t, x, v) , 
\ F{t, X, w) := x{t, x)d^Vl^m{t, x, w) + X{t, x, w)V^f{t, x, w) , 
Multiplying each side of the interaction identity (8) by 1r^ {v—w) and integrating 
both sides of the resulting equality in the variables v and w, 
(36) 

S:= [[ {AE-BD){t,x)dxdt 

= // ^R+iv — w) // C{s,x,v) I E{t,x,w)dt\ dxdsdvdw 
J J JJrxr \Js J 

+ 1r+(^' — w) // F(t,x,w) i / B{s,x,v)ds] dxdtdvdw . 

J J J JrxR \J-oo J 

As before, the right hand side of this identity involves terms of two different kinds: 
(37) 

Si = J J 1r^{v — w) J J X{s,x,v)V^ f{s,x,v) (^j E{t,x,w)dt^ dxdsdvdw , 

the other being 
(38) 

82= 1r+ {v — w) / / x(*) x)dv'D^m{s, x,v) I / E{t, x, w)dt I dxdsdvdw . 
J J J JnxR \Js J 

The inner integral 

/ E{t,x,w)dt= V^^{x{t,x)a'{w)f{t,x,w))dt- f{t + h,x,w)V'lx{t,x)dt 

J s J S J X 

s-\-h poo I \ 

X{t,y)a'{w)f{t,y,w)dy- a'{w)f{t + h,x,w){ d^x(.t + r,x)dT] dt , 



so that 



(39) 



/ E{t,x, 

J s 



w)dt 



< (llxlk~ + \\dtX\W)h\a' {w)\l[_u,u](w) ■ 
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Thus 

\Si\<Lih, 

with 

(llxlk- + Wdtxh^) J j \\X{-r,v)\\LAa'{w)\l^_u,u]{v)\^-u,u]i^)dvdw 
< (llxlk- + \\dtX\\L^)\\X\\i^A\a'U.(_u,u) =: • 

Moreover 

82= 'i-n+iv — 'w) // xis,x)dy'Dgm{s,x,v) I / E{t,x,w)dt\ dxdsdvdw 
J J JJrxr \Js J 

= — // x(t, x)I?Jm(s,x,«;) I / E{t,x,w)dt\ dxdsdw , 
JnJJnxR \Js ) 

so that 

\Si\<L2h 

with ^ ^ 

L2:=2(||x||L~ + ||aa||L0llxl|L- / V / |o'||m|. 

Collecting the bounds on Si and 52, one has 
(40) |5| <2(Li+L2)/i. 

On the other hand, according to Lemma 4.3 

(■T f 

x{t, x)'^A{u{t, x),u{t + h, x))dxdt 



(41) 



/O JR 

32 f-T 



{fi+m + 2)],, 7k 

Putting together (41) and (40) shows that 

which concludes the proof. □ 
Proof of Lemma 4-3. Assume that u < u. Then 

f l(fi,«](i^) if u > u > 0, 

Mu{v) - Mu{v) = I l[u,u]^v) + 1q{v) ifu>0>u, 
[ l[fi,„)(w) if > u > u. 

Therefore, \i — V < u < u <V , one has 

A{u,u) = jj l^^{v - w){a\v) - a\w))l(^u,u){v)l(u,u){w)dvdw 

{a'{v) — a{w))dv J dw 

> av / I {v — w)^dvdw 



U / fU 



u \Jw 

U PU 



'2 



(/3 + l)(/3 + 2) 



{u - uf+^ , 
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which is the announced lower bound when u < u. 
On the other hand, 

{M^{v) - Mu{v)){Mu{w) - Mu{w)) = {Mu{v) - Mu{v)){Mu{w) - M^iw)) 
for all u, u,v,w G R, so that 

A(u, u) = A(u, u) . 

With the previous inequality, this establishes the announced lower bound for all 

M, u G R. □ 

Observe that, for all uai G R 

{Mu{v) - Mu{v)){Mu{w) - Mu{w)) 

> l(inf(a,u),sup(a,«))('^)l(inf(fi,M),sup({t,M))(w') > 

for all f , w G R, so that Mu is an example of function satisfying the assumption 
(11) above. 

5. Regularizing Effect with One Convex Entropy 

In the case where the entropy condition is known to hold for one convex en- 
tropy, and with an entropy production rate that is a Radon measure with possibly 
undefinite sign, the following result was obtained by the first author: 

Theorem 5.1. [15, 16] Let a G C^(R) satisfy (10), and let ji he a signed Radon 
measure on R+ x R. Assume that the Cauchy problem for the scalar conservation 
law (25) has a weak solution u G L°°(R+ x R) satisfying 

(42) dtr,{u) + d^q{u) = -/x 

for some rj e (R) such that there exists /3' > 1 for which, given any ^ > 0, there 

exists r/ov > such that 

(43) v'i''^) ~ v'i^) ^ %,y(''^ ~ , whenever — V<w<v<V, 
with q defined by (4)- Then 

" e b^/;;„p^(r; X R) wtthp = p+p' + 2. 

More precisely, for all e > and each ^ G [— e, e] 

< (^M, + M2 + 2UMLA\\x\\L- + \\d,x\\L^)£ '^l'"') 
and, for each r G [0, e], 



< 



Nr + N2 + 2U\\x\\LA\\x\\L'>o + \\dtx\W) ( ^ I d\^l 

Jo J-R 



where U = ||x||l~, while Mi, M^, Ni, N2 are defined in (44), (45), (47) and (48) 
respectively and T,R> are chosen so that supp(x) C [0,T] x [—R,R]. 
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The proof of this rcsuh in [15, 16] is based on an argument that is reminiscent 
of Tartar's coinpcinsatcid compactness method for proving the convergence of the 
vanishing viscosity method for scalar conservation laws [39]. 

Whenever a,?7 G C^(R) satisfy 

a"{v) > a > and ri"{v) >r]o>0 for all u e R , 
one can take /3 = /3' = 1 in (10) and (43), and Theorem 5.1 predicts that u e 

In any case, the regularity obtained in Theorem 5.1 in the case where the weak 
solution u satisfies only one entropy relation with entropy production rate that is 
a signed Radon measure also belongs to the DcLellis-Westdickenberg [11] optimal 
regularity class i.e. the spaces bI^^^^^ for p > 3. 

A remarkable result due to Panov [31] states that, if m G L°°(R+ x R) is a weak 
solution of (25) with a e C^(R) such that a" > satisfying one entropy condition 
(42) with 1] G C^(R) such that 77" > and /i > 0, then it is the unique entropy 
solution of (25). In particular, it satisfies (42) for all 77 G C^(R) such that rj" > 0, 
with nonncgativc entropy production /i. Panov's result was subsequently somewhat 
generalized by DeLellis, Otto and Westdickenberg [10]. However, we do not know 
whether any weak solution of (25) with a G C^(R) such that a" > satisfying one 
entropy condition (42) with rj G C^(R) such that rj" > and fi a signed Radon 
measure must satisfy (42) for all convex entropies, i.e. whether such a solution 
satisfies the assumptions of Theorem 4.1. Therefore, Theorem 5.1 seems to be of 
independent interest. 

Below, we give a new proof of Theorem 5.1 based on the interaction identity of 
section 2 instead of the variant of compensated compactness used in [15, 16] . 

Lemma 5.2. Assuming that the functions a, t] and q belong to C^(R,) and satisfy 
(10) and (43) while q' = a'r]' , one has, for each V > 

{w - v){q{'w) - q{v)) - {a{w) - a{v)){r]{w) - r]{v)) 

ayTlovifi + 13') _ ,0+0' +2 

- (^ + /3' + l)(/3 + /3' + 2)' ' 

whenever v,w £ [—V, V] . 

In the most general case where a and 77 are convex functions, the inequality 

{w — v){q{w) — q{v)) — {a{w) — a{v)){r}{'w) — r]{v)) > , v,w gR. 

is stated without proof by L. Tartar: see Remark 30 in [39]. 
For a proof in the case where a, 77 G C^(R) with 

r]"{v) > ?7o > and a"{v) >a>0 for all w G R , 

corresponding with assumptions (10) and (43) with /3 = /3' = 1, see Lemma 2.3 in 
[15]. 
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Proof. One has 

{w — v){q{w) — q{v)) — {a{w) — a{v)){r]{w) — rj^v)) 



V J V 

W fW 



J 

J V J V 

(/3+/3'+l)(/3+/3'+2) I"' '^l 



□ 



Proof of Theorem 5.1. Pick x e C^(R+ x R), with support in ]0,T] x [-R,R], 

and let U := ||u||l-(r^ xR) • 

As above, we first estabhsh the regularity in the space variable x. Set 

A{t, x) -.^ x{t, x)V^u{t, x) , D{t, x) := x(i, x)V^r^{u){t, x) . 

Since u = f{t,x) satisfies (25)-(42), the function A and D defined above satisfy (6) 
with 



B{t,x) ■.= x{t,x)V^aiu){t,x), E{t,x) := x{t,x)V^q{u){t,x) , 



and 



r C{t, x) := V^u{t, x)dtx{t, x) + V^a{u){t, x)d^x{t, x) , 
\ F{t, x) := V^r]{u){t, x)dtx{t, x) + V^q{u){t, x)d^xit, x) - x(i, x)V^n . 
Using the identity (7) shows that 

X{t, zy{V^,uV^^q{u) - V^^a{u)V^^r^{u)){t, z)dzdt 



RxR 



= - 1 1 C{t,x) ( I D{t,y)dy] dxdt 

J J RxR \Jx ) 

- [f F{t, y) ( r A{t, x)dx] dydt . 
JJrxr \J-(x J 

Proceeding as in the proof of Theorem 4.1, we see that 

/y rv rv 

A{t,x)dx= V^{xu){t,x)dx - u{t,x + h)V^x{t,x)dx 
-OO J —OO J —QO 

ry+h ry f \ 

/ X^it, x)dx — / u{t, x + h)\ j dxX{tj X + z)dz 1 dx , 

Jy J -OO \Jo I 



so that 



Likewise 



/y 
A{t, x)d: 
-OO 

poo 

/ D{t,y)dy 

J X 



<INkoo(||xlk~ + l|5.xlU0N- 

< \H^)\\M\\x\\L<^ + \\d.x\W)\h\. 
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Therefore, assuming that \h\ < e, one has 

// C{t,x)(f D{t,y)dy] dxdt 
JJnxn \Jx J 



<Mi\h\, 



with 



(44) Ml := {2\\uU^\\dtx\W + l|aWI|L~l|5.x|Ui)h(«)||L~(||xl|L~ + . 
while 

<{M2 + M3)\h\, 



[[ F{t,y)( r A{t,x)dx)dydt 
JJuxh \J-oo / 



with 



(45) M2 := {2Uu)hoo\\dtx\\L^ + ||g(n)|Uoo ||a,xlUOII«IU~(llxlUoo + Wd^xh^) 
and 



Ms:=2\\uhA\\x\\L^ + \\d.x\W)\\x\\i 



T f-R+e 



J-R-e 



Thus 



/ / X{t, zf{V'iuVlq{u) - Via{u)V^ri{u)){t, z)dzdt 
J Jrxb. 

On the other hand, by Lemma 5.2 



< (Mi+M2 + M3)|/i| 



auVouil3 + P') 



(46) {V^^uV^^qiu) - V^Mu)T>'M^^)) > + ^, + ^^(^ + ^, + 2) 
so that the inequality above entails the estimate 

// x{t,zf\Viuf+^'+\t,z)dzdt 

JjRxR 

showing that 

-eLf+^'-^(R;;BV£r+^)'^+^'+^R)). 

As for the time regularity, pick h e [0, e] and set 

A{t, x) := x{t, x)V^u{t, x) , D{t, x) := x{t, x)V^r]{u){t, x) . 

Since u = f{t,x) satisfies (25)-(42), the function A and D defined above satisfy (6) 
with 

B{t, x) := x{t, x)V^a{u){t, x) , E{t, x) := x{t, x)V^q{u){t, x) , 

and 

r Cit, x) := V^luit, x)dtx{t, x) + V^a{u){t, x)d^x{t, x) , 

1 F{t, x) := V'lr]{u){t, x)dtx{t, x) + V'^q{u){t, x)d^x{t, x) - x{t, x)V'^ti ■ 
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At this point, we use tlie identity (8) which shows that 

/ / x{t, zf{V'iuVlq{u) - Vla{u)Vln{u)){t, z)dzdt 

= [[ C{s,x)( [ D{t,x)dt) dxds 
JJnxn \Js J 

+ 11 F(t,x)([ A{s,x)d.^ dxdt. 

J jRxR \J-oQ J 



As above 



f A{s,x)ds= I Vl{xu){s,x)ds- I u{s + h,x)V'lx{s,x)ds 



rt+h 



J xu{s:x)ds— I u{s + h,x)[ I dxx{s + T, x)dT \ ds , 



so that 



Likewise 



/A{s, x)d 
-oo 

D{t, x)dt 



<\\u\\LA\\x\\L^ + \\dtX\W)h. 

< II^WI|L~(||xl|L~ + ||5tX||L~)/l. 



Therefore 



RxR 



C{s,x) ( / D{t,x)ds dxdt 



<N^\h\ 



with 



(47) iVi := {2\\uUAdtx\W + llaWIU'^ P^xlUOhMlk-dlxlk- - 
while 

ft' / /■* \ 

< {M2+M^)\h\. 



1 1 F{t, x)( I A{s, x)ds \ dxdt 



with 



(48) 7V2 := (2||r?(M)||ioc||atx||Li + l|9(«)lk~ l|a.xlUOII«lk~ (llxlkoo + Wd^xU^) 



and 



N^ := 2||«|U»(||x||l~ + H^txIUOIIxlk- T^' t N 

>/0 J-R 

Using the inequaUty (46), we conclude that 

x{t,zf\D\uf^f^'^\t,z)dzdt 



< 



'RxR 

{l3 + (3' + l){P + l3' + 2) 
auVou{P + P') 



{Ni+N2 + N3)h, 



so that 



,^+/3'+2/„. pl/(/3+/3'+2),/3+/3'+2.„* ^^ 
^ S hoc l**-' ^oo,loc 

which concludes the proof. 



□ 



OPTIMAL REGULARITY FOR CONSERVATION LAWS 
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